CLNS 02/1788 
CLEO 02-8 

Measurement of B(B~— >D 7r _ ) and B(B ^D + tt~) and Isospin 

Analysis of B^Dtt Decays 

CLEO Collaboration 
(Dated: February 4, 2008) 

Abstract 

We present new measurements of branching fractions for the color- favored decays B~^D°tt~ 

O ■ _ 

and Tr^D+it-. Using 9.67 x 10 6 BB pairs collected with the CLEO detector, we obtain the 



branching fractions B(B-^D°tt-) = (49.7 ± 1.2 ± 2.9 ± 2.2) x 10" 4 and B{B^ — >D^tt~) = (26.8 ± 
1.2 ± 2.4 ± 1.2) x 10~ 4 . The first error is statistical, the second is systematic, and the third is due 
to the experimental uncertainty on the production ratio of charged and neutral B mesons in T(AS) 
decays. These results, together with the current world average for the color-suppressed branching 
O ■ fraction B(B^^D°tt°), are used to determine the cosine of the strong phase difference Si between 
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the I = 1/2 and I = 3/2 isospin amplitudes. We find cos 5/ = 0.863I^23Tao35-ao30> and obtain a 
90% confidence interval of 16.5° < 5j < 38.1°. This non-zero value of 5i suggests the presence of 
final state interactions in the Dir system. 
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This paper presents the results of measurements of the branching fractions for B~^D n~ 
and B°-^D + -it~ and the extraction of the strong phase difference 5i between the 1 = 1/2 and 
1 = 3/2 isospin amplitudes in the Dir system. These decays are an excellent testing ground 
for the theoretical description of hadronic 5-meson decays. Our understanding of these 
decays has improved considerably during the past few years with the development and ap- 
plication of Heavy Quark Effective Theory (HQET) [[l], ||| and Soft Collinear Effective Theory 
(SCET) 0. Originating from the simple, but very effective, idea of color-transparency [|J, 
the factorization hypothesis has been put on a more solid basis, and in the case of B^Dir, 
has been proven within the framework of SCET. 

The recent observation |5|, || of the color-suppressed 5°— >D°ti decay 1 completed the 
measurement of the Dir final states and was used to determine the cosine of the strong 
phase difference cos 5/ = 0.89 ± 0.08, a value which is consistent with one. A value of cos Si 
inconsistent with one would signal the presence of final-state interactions in the B^Dn 
process 0, §f. In this paper, we present improved measurements of branching fractions for 
the color-favored decays B~^D°tt~ and B°— >D + tt~ based on a larger data set than that 
from which the previous results were obtained, as well as a new evaluation of cos 5j which 
takes into account the correlations among the various contributions to the overall systematic 
error. 

This analysis uses e + e~ annihilation data recorded with the CLEO detector at the Cornell 
Electron Storage Ring. The integrated luminosity of the data sample is 9.15 fb^ 1 collected 
on the T(AS) (on-resonance), corresponding to 9.67 x 10 6 BB pairs, and 4.35 fb^ 1 collected 
60 MeV below the BB threshold (off-resonance), which is used for background studies. The 
results we present in this paper for B(B~^D -k~) and B(B°^D + n^) supersede those in 
the CLEO publication, Ref. 0, which were based on a 1.3 fb^ 1 subset of the data used in 



the present analysis. Data were recorded with two detector configurations, CLEO II 1 1 



and CLEO II. V Cylindrical drift chambers in a 1.5 T solenoidal magnetic field measure 
momentum and specific ionization (dE/dx) of charged particles. Photons are detected using 
a CsI(Tl) crystal electromagnetic calorimeter, consisting of a barrel-shaped central part of 
6144 crystals and 1656 crystals in the forward regions of the detector (endcaps). In the 



1 Throughout this paper, charge conjugation is implied. 
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CLEO II. V configuration, the innermost tracking chamber was replaced by a three-layer, 
double-sided silicon microvertex detector, and the main drift chamber gas was changed from 
argon-ethane to a helium-propane mixture. 

In our analysis, we impose quality requirements on charged particle tracks and improve 
the purity of pion and kaon used to reconstruct D mesons by using dE/dx information if 
the particle momentum is less than 800 MeV/c. The neutral D mesons are reconstructed 
using three decay modes: K~ir + , K~ir + Tr° and K~tt + 7t~tt + . Charged D mesons are similarly 
reconstructed via the mode K~tt + tt + . In each case, D meson candidates are required to have 



a mass within 3cr (standard deviations) of the PDG D mass |J2] before kinematic fitting. 
Resolutions for the various D modes range from 6 to 12 MeV. 

Each B meson candidate is reconstructed using the four-momentum of the mass- 
constrained D meson and an additional charged track in the event (assumed to be a pion). 



Candidates are then identified using the beam-constrained mass M B = J E^ eam — p 2 B , where 
-Sbeam denotes the beam energy and ps the candidate momentum, and the energy differ- 
ence AE defined by AE = Ed + E n — E^ cam , where Ed and E n are the D meson and tt 
energies, respectively. Preselection of B candidates requires Mb > 5.24 GeV/c 2 and AE 



to be between —50 and 50 MeV. Additionally, we calculate the sphericity vectors [13|] of 
the B daughter particles and of the rest of the event. We require the absolute value of the 
cosine of the angle between these two vectors to be less than 0.8. The distribution of this 
angle is strongly peaked at ±1 for continuum background and is nearly flat for BB events. 
We also require events to satisfy R2 < 0.45, where R2 is the ratio of the second to zeroth 



Fox- Wolfram moments of the event ||14|| . Finally, for events with more than one B meson 
candidate, the candidate with the smallest |Ai?| is chosen. 

To obtain event yields for B^Dn~ for each D meson decay mode, the Mb distribution 
of candidates surviving the above slection cuts are fit using a binned maximum likelihood 
fit. The function used is a Gaussian for the signal plus an empirical background function, 



f(M B ) = AM BV /l-(M B /E bcam ) 2 expa(l-(M B / J E; bcam ) 2 ), having a fixed E hcam = 5.29 
GeV. All other parameters in both background and signal functions are allowed to float in 
the fit. The fitted Mb distributions for each of the D meson decay modes are presented in 
Fig. § 

A small, non-negligible background from the decay B^DK~ contributes to the yields 
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FIG. 1: Fitted M B distributions in for a) Br ^D q {K--k + )it- b) B~ ^D°{K-ir + ir Q )ir c) 
B-^D (K--K + TT--K + )7r- and d)U ^D+(if-7r+7r+)7r-. 

obtained by the fit procedure described above. We have, therefore, simulated this back- 
ground via Monte Carlo to determine the fraction of feed-through to the Dtt~ sam- 



ple, and performed a subtraction using the average of the two measurements [15, 16 



of B( y B-^D K-)/B(B-^D°n-) = 0.071 ± 0.009 and the recent measurement (T§ of 



B(B°^D+K-)/B(B°^D + it-) = 0.068 ± 0.017. The amount of DK feed-through is found 
to be approximately (4 ± 1)% of the Dtt yield. We then reduce the event yields obtained in 
the fit to the data by this fraction. 

Using efficiencies determined by applying the above method of analysis to samples of signal 
Monte Carlo events, we obtain the branching fractions for the processes under investigation 
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TABLE I: Results for the branching fractions B(B-^D°tt-) and B(lf^D + ir-). Fit yields with 
errors and efficiencies are obtained as described in the text. The errors given for the efficiencies 
correspond to the Monte Carlo statistical errors for each mode. The D mode branching fractions 



and the branching fraction errors have been taken from the PDG [12]. The errors reported for the 
measured B branching fraction are the statistical errors only. The current PDG average values for 
the two branching fractions have been included for comparison. 



B-^D°ir- 


D° Decay Mode 


Yield Efficiency (%) D° mode B(%) 


B(B-^D°tt-)(x10- 3 ) 


K-TT+ 
K-1T+1T 

K~ir + ir~7r + 


820 ±31 45.4 ±0.3 3.83 ± 0.09 
1200 ±45 17.1 ±0.2 13.9 ±0.9 
740 ± 33 20.9 ± 0.3 7.49 ± 0.31 


4.90 ±0.18 
5.20 ±0.19 

4.91 ±0.22 


PDG 




5.3 ±0.5 


B°^D+tt- 


D + Decay Mode 


Yield Efficiency (%) D+ mode B(%) 


B(B°^D+ir-)(xW- 3 ) 




764 ±33 32.8 ±0.4 9.0 ± 0.6 


2.68 ±0.12 


PDG 




3.0 ±0.4 



from the event yields corrected for the DK feed-through: 

Corrected Yield 



B(B^Dtt) 



(1) 



e x B(D^f.s.) x N(T(AS)) x 2/' 
where / represents /+_ or / o, the charged or neutral B meson production ratios at the T(4S r ), 
as appropriate. The corrected yields, efficiencies and final branching fraction obtained for 
each D decay mode are shown in Table |. We have assumed / + _ = /oo = 0.5. 

The three B^Dir decay branching fractions (the two color-favored modes, which we re- 
port measurements for in the present paper, as well as the color-suppressed mo de B —>D° 
form a complete set of branching fractions with which we may calculate cos Si, the cosine of 
the strong phase angle difference between the two isospin amplitudes 1=1/2 and 7 = 3/2 
that contribute to the decay process. The expression for cos 5/, following Ref. |7|], is: 

3T(D-n + ) + T(D\ + ) - QT(D°ir°) 



COS 5; 



4|^l/ 2 ^3/ 2 | 
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(2) 



where the isospin amplitudes A3/2 and A1/2 are given by 



\A 3/2 \ 2 = T(D\ + ), and (3) 

IA /2 | 2 = ~ (rp-vr + ) + r(DV)) -§r(EV). (4) 



The calculation of cos Si in the Dir system takes into account correlations of system- 
atic errors between the two color-favored decay modes B~^D°tt~ and B°-^D + ti~ . It also 
considers the fact that some of the systematic errors in the measurement of B(B~^D°tt~) 
using the three D° decay modes are correlated. Further, apart from the errors on /oo and 

/_i (which are anticorrelated) , we treat the errors between the two color-favored B^Dit 

modes and the color-suppressed fi°^D 7r° mode as uncorrelated. This treatment is justified 
since the systematic error on the color-suppressed mode is dominated by the background 
parameterization and fit uncertainties, whereas such contributions are not dominant for the 
color-favored modes J5|, |6|. 

We estimate the following systematic error contributions to our results for these mea- 
surements: 1% per track for track finding and fitting, 2% for the total number of BB pairs, 
2% per track for which dE/dx is used, 2.5% for the cuts used in the analysis and 1% for 
the DK feed-through subtraction. Other systematic errors include 2% for 7r° finding in the 
case of the D—>-K~ir + 7r submode, 2.3% — 7% for D-meson branching fractions, 2 — 3% for 
background parameterization and fitting, and 0.7% — 2% for Monte Carlo statistics. The 



experimental errors of 4.5% on the individual quantities / + _ and / o HH are reported as a 
separate systematic error in our final result. 

The overall systematic error for our measurement of B(B° ^D + n~) is obtained by stan- 
dard error propagation of the individual contributions. However, in order to extract the 
correct overall systematic errors for B(B~^D°tt~) and for cos^j-, we must take into account 
the correlation among the systematic error contributions for each of the D submodes. To do 
this, we perform Monte Carlo experiments in which we vary the measured branching frac- 
tions by their various systematic errors. In each experiment and for each systematic error 
contribution, we generate multiplicative correction factors according to a Gaussian distri- 
bution. The combined B~-^D°tt~ branching fraction and cos 5/ are then calculated from 
the values which have been varied as described above for each Monte Carlo experiment. 
From the complete ensemble of 2 x 10 6 Monte Carlo experiments, we obtain the probability 



distribution functions and errors for B(B ^D°tt ) and cos 5/, which are shown in Fig. |^. 
We thus obtain the following final results for the color-favored branching fractions. 

B(B-^D ir~) = (49.7 ± 1.2 ± 2.9 ± 2.2) x 10 -4 , 
B(B°^D + n-) = (26.8 ± 1.2 ±2.4 ±1.2) x 10~ 4 

In each measured quantity, the first error is statistical and the second is systematic. The 
third error is a separate systematic error which corresponds to the experimental uncertainty 
of the production fraction of charged (or neutral, as appropriate) B mesons in T(4S) decays. 

Our results for B(B°^D + 7t~) and for B(B~—>-D 'n-~) each reflect improvement with re- 
spect to the present PDG average values |12|]. Our result for £>(£>°— >/} + 7r~) may be directly 
compared with the prediction of Ref. |1|] for this decay. Their prediction of 32.7 x 10~ 4 is 
marginally consistent with our result. 

The largest contribution to the overall systematic error in our result for B(B°-^D + n^) 
is the 6.6% relative systematic error due to the D branching fraction. We therefore re- 
port the following result which is independent of the D + — >K~7t + tt + branching fraction: 
B(B°^D + 7i-) x B(D + ^K-n + n + ) = (2.41 ± 0.11 ± 0.15 ± 0.11) x 10~ 4 . 

Using the CLEO measurement || of the color-suppressed branching fraction, 
B(B°^D°it ) = 2.74l{j;32 ± 0.55 x 10" 4 , and the PDG(2002) ratio of B lifetimes, 
t(B + )/t(B°) = 1.083 ± 0.017, we obtain 

cos = 0.877 ±0.0301^1^3?- 

The error distributions derived from the ensemble of Monte Carlo experiments for cos Si 
and 5 1 are shown in Fig. ^|. Integrating the 5i distribution over the physical region | cos5/| < 
1, we obtain a 90% confidence interval: 

13.6° <5 r < 38.3°. 

Our final results for cos 5/ and 5j are based on the average of both measurements of 
B(B°^D°n ) = 2.92 ± 0.45 x 10~ 4 § §. Using this average, we obtain 

m „l _ n Q«q+0. 024+0.036+0.038 
LVbUi — U.OUO_ .023-0 .035-0.030" 

Similarly, we obtain our final result for Si, a 90% confidence interval of 

16.5° < 5/ < 38.1°. 
9 



Using our results for B(B —^D°tt ) and B(B ^D + n ), we also calculate the ratio of the 
I = 1/2 and J = 3/2 isospin amplitudes, A1/2/A3/2 = 0.69 ± 0.03 ± 0.06 ± 0.06. In the heavy 
quark limit, A1/2/A3/2 = l- 2 Corrections to this are C(Aqcd/^c), which is consistent with 
our result. 

In summary, we have measured the branching ratios for the color-favored B—>Dtt de- 
cays, and used these measurements, together with the current average of measurements of 
B(B ^D°n ), to determine the value of the cosine of the strong phase difference 5i in the Dir 
system, and the ratio of I = 3/2 and 1 = 1/2 isospin amplitudes. Our result for cos Si differs 
from one by approximately 2.3cx and thus suggests the presence of final-state interactions in 
B^Dn decays. 

We thank Jonathan Rosner and Matthias Neubert for useful discussions. We gratefully 
acknowledge the efforts of the CESR staff in providing us with excellent luminosity and 
running conditions. M. Selen thanks the PFF program of the NSF and the Research Corpo- 
ration, and A.H. Mahmood thanks the Texas Advanced Research Program. This work was 
supported by the National Science Foundation, and the U.S. Department of Energy. 



2 The ratio given here is based on the formalism of Ref.0. It is equivalent to A1/2/ {V2A3/2) according to 
that of Ref. D 
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FIG. 2: The error distributions for cos 5i and 5i obtained from the ensemble of 2 x 10 Monte Carlo 
experiments described in the text. The shaded area in the cos 6i plots is the ±lcr window (the 90% 
C.L. region in the 5i plots). The upper two plots show the distributions for a) cos 5/ and b) 5i 
obtained using only the CLEO measurement of B(B°^D ir ). The lower two plots are distributions 
for c) cos 5[ and d) 5[ obtained using both CLEO and Belle measurements of B(B°^D ir°). 
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